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We present binary galaxy merger simulations of gas-rich disks (Sp-Sp), of early-type galaxies and disks (E-Sp, mixed 
mergers), and mergers of early-type galaxies (E-E, dry mergers) with varying mass ratios and different progenitor mor- 
phologies. The simulations include radiative cooling, star formation and black hole (BH) accretion and the associated 
feedback processes. We find for Sp-Sp mergers, that the peak star formation rate and BH accretion rate decrease and the 
growth timescales of the central black holes and newly formed stars increase with higher progenitor mass ratios. The ter- 
mination of star formation by BH feedback in disk mergers is significantly less important for higher progenitor mass ratios 
(e.g. 3:1 and higher). In addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers 
and mixed E-Sp mergers. 
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1 Introduction 

Recent large-scale statistical galaxy surveys (e.g. the Sloan 
Digital Sky Survey (SDSS) and the two-degree Field (2df) 
survey) have revealed a robust bimodality in the observed 
galaxy population. Galaxies above a critical stellar mass of 
M cr it ^ 3 x 1O 1O M are typically non-star forming red 
spheroidal systems with old stellar populations that predom- 
inately live in dense environments, whereas galaxies below 
this critical mass are typically blue, star-forming disk galax- 
ies that lie in the field (e.g. Bell et al. 2003; Kauffmann et 
al. 2003; Baldry et al. 2004; Balogh et al. 2004). 

The observed bimodality can also be seen as a mani- 
festation of cosmic downsizing, in which the most massive 
galaxies formed a significant proportion of their stars at high 
redshifts above z > 2. The star formation is efficiently 
quenched in these systems by z = 1 as manifested in the 
observed population of extremely red non-starforming mas- 
sive ellipticals (EROs, e.g. McCarthy 2004; Vaisanen & Jo- 
hansson 2004). The lower mass systems, on the contrary, 
have typically been forming stars throughout the whole cos- 
mic epoch (e.g. Glazebrook et al. 2004; Juneau et al. 2005; 
Cimatti, Daddi, Renzini 2006). 

The quenching mechanism responsible for the observed 
termination of star formation in massive galaxies at redshifts 
below z < 2 — 3 needs to be both energetic enough to trig- 
ger the quenching and long-lasting enough to maintain the 
quenching over a Hubble time. There exists several theoret- 
ical explanations for the quenching, including the feedback 
from AGNs (e.g. Bower et al. 2006; Croton et al. 2006), 
gaseous major mergers triggering star burst and/or quasar 
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activity (e.g Naab, Jesseit, Burkert, 2006), quenching by 
shockheated gas above a critical halo mass (Dekel & Birn- 
boim 2006; Birnboim, Dekel, Neistein 2007) and gravita- 
tional quenching by clumpy accretion (Naab et al. 2007; 
Dekel & Birnboim 2008). 

In this paper we study the termination of star forma- 
tion by black hole (BH) feedback in binary galaxy merg- 
ers. The study by Springel, Di Matteo, Hernquist (2005a) 
showed that BH feedback efficiently terminates star forma- 
tion in equal-mass disk mergers. Here we confirm this result 
and expand the analysis to include unequal-mass disk merg- 
ers, mixed E-Sp mergers and dry E-E mergers. This paper is 
structured as follows. In ^2] we briefly describe our simula- 
tion code and galaxy merger setup. We discuss the effect of 
varying the initial gas fractions, orbits, merger mass ratios 
and galaxy progenitors for terminating star formation in ^3] 
Finally, we summarize and discuss our findings ^4j 



2 Simulations 

We perform the simulations using the parallel TreeSPH- 
code GADGET-2 (Springel 2005) including star formation 
and the associated supernova feedback, where the multi- 
phase interstellar medium (Efstathiou 2000; Johansson & 
Efstathiou 2006) is modeled using the sub-resolution mo- 
del developed by Springel & Hernquist (2003). We imple- 
mented a black hole feedback description similar to the Sp- 
ringel et al. (2005b) model into our version of the code. In 
this model the BH accretion rate is parameterized using a 
Bondi-Hoyle-Lyttleton description, with the maximum ac- 
cretion rate limited to the Eddington rate. Finally, some frac- 
tion of the accreted rest mass energy is available as thermal 
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Fig. 1 The total star formation rate (top), the total black 
hole accretion rate (middle) and the evolution of the total 
black hole mass (bottom) as a function of time for three 3:1 
(solid lines) and one 1 : 1 (dashed line) merger with initial 
gas mass fractions of 20% (black), 40% (green) and 80% 
(red) (left panel) and for three 3:1 mergers with varying ini- 
tial gas mass fractions (black, green, red lines) and three 
3:1 mergers with varying orbital and initial geometries for a 
fixed gas fraction (blue lines) ( right panel). The filled circles 
indicate the time of merging of the BHs. 

feedback energy that couples to the surrounding gas. Fur- 
ther details about the code implementation and parameter 
choices can be found in Johansson, Naab, Burkert (2008), 
J08 hereafter. 

The disk galaxies are setup as described in J08, with 
Hernquist dark matter profiles populated by exponential dis- 
ks with initial gas fractions ranging from 20% to 80%. All 
models contain initially a bulge with a third of the total disk 
mass. We also model mergers of early-type galaxies, where 
the early-type galaxies are setup initially using merger rem- 
nants of binary disk mergers (Naab & Burkert 2003). We 
set the primary galaxy models up with 20,000 gas and stellar 
disk particles, 10,000 bulge particles and 30,000 dark matter 
particles and scale the other models correspondingly. The 
gravitational softening lengths of gas, newly formed stars 
and the black hole particles are set to e = 0.1/i _1 kpc and 
the softening lengths of the disk, bulge and more massive 
dark matter particles are scaled with the square root of the 
particle masses resulting in e = 0.2/i _1 kpc for the bulge 
and disk particles and e = 0.8/i _1 kpc for the dark matter 
particles, respectively. 

The galaxies are initially placed on parabolic orbits wh- 
ere the initial separation of the progenitors is set to the mean 
of the two galaxy virial radii and the pericentric distance to 
the mean of the two disk scale radii. All simulations pre- 
sented in this paper were evolved for a total of t = 3 Gyr 



using the local Altix 3700 Bx2 machine hosted at the Uni- 
versity Observatory in Munich. 

3 Results 

In the left panel of Fig. \T\ we plot the total star formation 
rates, BH accretion rates and the BH mass growth for three 
3:1 mergers (solid lines) with 20% (black), 40% (green), 
80% (red) initial gas mass fractions together with one 20% 
gas fraction 1 : 1 merger (dashed lines). The star formation is 
very efficiently terminated by the BH feedback in 1 : 1 merg- 
ers, compared to a generally much shallower decline in star 
formation for 3:1 mergers. In addition, the final BH mass 
are lower in 3:1 mergers typically by a factor of 2-5, but 
with a relatively large scatter depending on the progenitor 
masses and initial gas fractions. 

In the right panel of Fig. [T] we show the results for three 
3:1 mergers with varying initial gas mass fractions for a 
fixed orbit and initial disk geometry (black, green, red lines) 
and for three 3:1 mergers with varying initial orbits and ori- 
entations for a fixed gas mass fraction (blue lines). The ini- 
tial gas fraction has a large effect on the height of the star 
formation and BH accretion peaks, with larger initial gas 
fractions producing higher values, as expected. This results 
also in relatively large differences in the final BH masses, 
with the /g as = 0.8 simulations producing final BH masses 
that are larger by an order of magnitude compared to the 
/gas = 0.2 runs. The variation of the orbit and initial ge- 
ometry for a fixed gas mass fraction produces much smaller 
differences. The peaks of the star formation rates and BH 
accretion rates only vary within a factor of two with chang- 
ing orbits and initial disk geometries. 

In Fig J2] we study the star formation and BH accretion 
histories for unequal-mass mergers with varying mass ra- 
tios. To this end we ran five mergers with mass ratios of 
1:1,2:1,3:1,4:1 and 6: 1 on co-planar prograde orbits with 
initial gas mass fractions of 20%. As can be seen in Fig. [2] 
increasing the mass ratio of the merger systematically low- 
ers the peak star formation rate and increases the duration 
of star formation activity after the merger. For the highest 
mass ratio merger the star formation rate is virtually con- 
stant throughout the simulation with only a mild peak dur- 
ing the first passage (see also di Matteo et al. 2007). The fi- 
nal BH masses are systematically lower for increasing mass 
ra-tio of the merger. Furthermore, the slope of the Mbh 
growth as a function of time becomes shallower with in- 
creasing merger mass ratio. There is also a systematic delay 
in the time of the BH merger with increasing mass ratio, as 
indicated by the filled circles in Fig. [2] For the lower mass 
ratio mergers the peak of the BH activity (the filled triangles 
in Fig. typically occurs shortly after the merging of the 
BHs, whereas for the higher mass ratio mergers the peak of 
the BH activity is not directly related to the merging time of 
the BHs. 

In Fig. [3] we study the duration of the BH accretion and 
star formation activity as a function of merger mass ratio. 
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Fig. 2 The evolution of the star formation rates, black 
hole accretion rates and BH mass as a function of time for 
co-planar prograde 1:1 (black), 2:1 (green), 3:1 (red), 4:1 
(blue) and 6:1 (orange) mergers. The filled circles indicate 
the time of the BH merger, the filled diamonds in the top 
panel and the filled triangles in the bottom two panels show 
the location of the maximum star formation and BH accre- 
tion peaks, respectively. 

We define a half-mass growth time AT x / 2 during which half 
of the final BH mass and half of the total stellar mass is 
formed respectively. In both cases A7\ / 2 is centered on the 
peak of the corresponding activity, the BH accretion on the 
maximum of the BH accretion rate (the triangles in Fig. [2]) 
and the star formation rate on the peak of the SFR marked 
with diamonds in Fig. [2 The resulting half-mass timescales 
are shown in the bottom panel of Fig. [51 The AT!/ 2bh is 
strongly correlated with the mass ratio of the merger. For 
1:1 and 2:1 mergers the growth of the Mbh is very con- 
centrated in time, with half of the final BH mass growth 
occurring in less than 100 Myr. For the 3:1 and 4:1 merg- 
ers AT!/2 BH ~ 0.5 Gyr, with the 6:1 merger resulting in 
AT 1/2bh ~ 1 Gyr. 

The corresponding stellar half-mass timescales ATi 

_ 2 star 

(triangles in the bottom panel Fig. [3]) also show a clear cor- 
relation with the mass ratio of the merger. In the 1 : 1 merger 
half of the final stellar mass is formed in a short major burst 
lasting about ATi ~ 0.15 Gyr. This value is compa- 
rable to the star formation timescale derived by Cox et al. 
(2008), who calculated a full width at half maximum of 
FWHM ~ 0.1 Gyr for the star formation peak of a typical 
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Fig. 3 The maximum BH accretion and star formation 
rates and timescales as a function of merger mass ratio. 
The bottom panel shows the half mass growth time of the 
BHs (circles) and the stellar mass (triangles) centered on 
the maximum BH accretion and star formation rates, respec- 
tively. The top panel shows the corresponding peak BH ac- 
cretion (circles) and star formation rates (triangles) during 
the corresponding half mass growth time A7\ / 2 . 

1 : 1 merger. For higher mass-ratio mergers the resulting star 
formation timescales are ATi ~ 0.7 — 1.0 Gyr, with 

2 star 

the highest mass-ratio merger having the longest timescale 
ofATi -1.5 Gyr. 

2 star 

In the top panel of Fig. [3] we plot the corresponding 
maximum black hole accretion and star formation rates. By 
defining the variable q as the mass ratio between the pri- 
mary and secondary component we can fit the logarithms 
of the peak BH accretion and star formation rates with the 
following linear relation 

log Maximum [MQyr -1 ] = ao + ai x q, (1) 

where ao and a\ are the inferred normalization and slope re- 
spectively. Both the maximum BH accretion rates and star 
formation rates are well fitted by Eq. [T] (dotted lines in Fig. 
[3]) resulting in the following best fitting parameters respec- 
tively: (a ,BHAR = 0.31,ai iB HAR = -0.47) (a ,sFR = 
1.82,ai y sFR = —0.24). The ratio of the peak BH accre- 
tion rate to the peak star formation rate is thus of the or- 
der of MBH,peak/M S F, P eak ~ 10~ 2 . On the other hand, 
the ratio of the mean BH accretion rate to the mean star 
formation averaged over the whole simulation is closer to 
M BH ,mean /M SF .mean ~ 10 3 , with this ratio varying sys- 
tematically between M B H,mean/M S F,mean = 2 • 10" 3 (1:1 

merger) and M B h, mean /M SF , mean = 0.2 • 10" 3 (6:1 mer- 
ger). Interestingly, recent observations of z ~ 2 galaxies by 
Daddi et al. 2007 also find indications of an universal ratio 
of Mbh/Msf ~ 10 -3 and this is also expected from the 
observed Mbh — M* relation. 
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Finally, we analyse in Fig. [4] the star formation and BH 
accretion rates for a sample of four mixed E-Sp mergers 
(left panel) and four dry E-E mergers (right panel). The star 
formation rate of the mixed E-Sp mergers is lower than in 
Sp-Sp (Fig. [T} mergers, due to the lower amount of cold 
gas available for star formation. The disk progenitor con- 
tains / gas = 20% of gas initially, whereas the early-type 
progenitors typically have an initial gas mass fraction of 
/gas ~ 5%, with typically only ~ 1% of this gas being 
cold and dense. After the merger of the E-Sp galaxies the 
star formation rate declines rapidly in all the merger rem- 
nants. We define a BH mass growth factor as /BH,insitu = 
AM B H,insitu/M B H,finai, which gives the ratio of BH mass 
growth due to gas accretion during the simulation with re- 
spect to the final BH mass. Quantitatively, the fraction of 
the BH mass that accumulates by gas accretion during the 
mixed mergers is in the range of /BH,insitu ~ 20 — 50%, 
with typical mean values of /BH,insitu ~ 30%. For the E- 
E mergers the initial star formation rates are generally very 
low due to the low gas fractions of / gas ~ 1 — 5%. The ini- 
tial gas fraction directly depends on the strength of the ini- 
tial interaction that gave rise to the merger remnants used as 
progenitors for the E-E remergers. Increasing the masses of 
the progenitors and using more direct planar orbits produces 
more violent initial encounters, thus decreasing / gas . The 
star formation is very effectively terminated shortly after the 
merging of the E-E progenitors on comparable timescales to 
1:1 Sp-Sp mergers and thus more efficiently than in 3:1 Sp- 
Sp and mixed E-Sp mergers. 

4 Conclusions 

In this paper we have studied the termination of star for- 
mation in merger simulations including BH feedback. We 
find that the termination of star formation by BH feedback 
is significantly less important for unequal-mass disk merg- 
ers compared to equal-mass disk mergers. The timescale 
for star formation termination systematically increases with 
increasing progenitor mass ratios. Similarly, a systematic 
increase is seen in the half-mass growth timescales of the 
BHs, with this timescale varying from ~ 0.1 Gyr for equal- 
mass mergers to ~ 1 Gyr for 6:1 mergers. This systematic 
trend can be used as input in modeling BH accretion more 
realistically in semi-analytic galaxy formation models (e.g. 
Croton et al. 2006). For mass-ratios of 3: 1 and higher merg- 
ers with BH feedback are unable to completely quench the 
star formation, with the merger remnants showing star for- 
mation rates roughly on the pre-merger level even 1 Gyr af- 
ter completion of the merger. In addition, the star formation 
is efficiently terminated in mixed E-Sp and dry E-E mergers 
due to the presence of the fully grown super-massive BHs 
in the early-type progenitors. 

Acknowledgements. The numerical simulations were performed 
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Cluster of Excellence: "Origin and Structure of the Universe". 
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Fig. 4 The total star formation rate (top), the total black 
hole accretion rate (middle) and the evolution of the total 
black hole mass (bottom) as a function of time for four E- 
Sp mixed mergers (left panel) and four E-E remergers ( right 
panel). The gas disks have initially / gas = 20% and the 
virial velocities as indicated in the Figure. The filled circles 
indicate the time of merging of the BHs. 
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